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Lepton flavour violation is considered in two models, Standard Model with additional heavy Dirac neutrinos 
and minimal supersymmetric standard model (MSSM) with 5*0(10) theoretical frame and minimal SUGRA 
(mSUGRA) SUSY-breaking mediation mechanism. The models are briefly explained. The structure of the 
amplitudes for neutrinoless decays of the lepton into a lepton and semileptonic neutrinoless decays of a lepton is 
given. The method for identification of meson content in the quark currents is exposed. The comparison of the 
amplitudes is made showing which of the processes have the largest branching ratios. The numerical results for 
the dominant decay rates are presented. 



1. Introduction 

Neutrino oscillation data not only show that 
the neutrinos do have mass and do mix, but 
these mixings and differences of the squares of 
the neutrino masses are well determined. Exten- 
sions of the Standard Model (SM) are therefore 
necessary, and any such extension must be able 
to describe the neutrino data together with all 
other experimental data which SM describes well. 
Two of such models are considered here. The 
simpler of these models, SM extended by heavy 
(quasi-)Dirac neutrinos '1"2'3' (called model A 
here) has only additional heavy (quasi-) Dirac 
neutrino fields compared to the SM field con- 
tent. The second model '4"5' 6l7l8ig] (model B) is 
based on the minimal supersymmetric standard 
model (MSSM) with heavy right-handed neutri- 
nos |1UI11| . with model parameters determined by 
the underlying GUT model and through the min- 
imal SUGRA (mSUGRA) SUSY-breaking medi- 
ation mechanism (12) . 

The existence of the neutrino mixings implies 
that the lepton flavours are violated. An indepen- 
dent confirmation of the lepton flavour violation 
(LFV) is necessary. The last improvements of the 



LFV upper bounds strongly motivate further 
theoretical investigation of these processes. The 
structure of the LFV amplitudes and a numerical 
study of the corresponding decay rates in model A 
and model B are the subject of of this paper. The 
main topic of this paper are semileptonic (SL) 
neutrinoless LFV decays of charged leptons. The 
LFV decays of the charged lepton into lepton and 
photon {£ — > £''-f) are considered also, to deter- 
mine the bounds on the parameters in the model 
B. 



2. Models 

The model A is inspired by a £^8 x Eg heterotic- 
superstring theory. Along with the SM neu- 
trino fields, it contains a number of right-handed 
neutrino fields and equal number of left-handed 
E'e-singlet neutrino fields Ill2ll4j . The right- 
handed neutrino fields mix through the Dirac- 
type Yukawa couplings with SM neutrino fields 
and right-handed neutrino fields. The mass ma- 
trix has three zero eigenvalues while the other 
eigenvalues are are large (typically ~ 1 GeV). 
The masses and the mixings of the light neutrinos 
may be induced through the small Majorana-type 
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Yukawa couplings of the left-handed Eg singlet 
neutrino fields 14 15 . Fitting of the neutrino 
oscillation data is easily achievable. The small 
Majorana mass terms promote the heavy Dirac 
fields to heavy quasi-Dirac fields. The structure 
of the mass matrix without Majorana terms as- 
sures conservation of the total lepton number jl4l 
llfillTj , but the individual lepton fiavours are vio- 
lated through the neutrino mixings |16ll7ll8ll9j . 
Through the neutrino mixings, the SM neutrino 
fields become massless (light) neutrino fields {i^i, 
i — 1,2,3) with small admixture of the heavy 
(quasi-)Dirac fields {Ni, ? = 4, • • •). Through the 
loops, in which only the massive neutrinos give 
a contribution, the heavy (quasi-)Dirac neutrino 
components of the SM neutrino field give rise to 
the LFV effects that are not GIM suppressed. 
The CKM-type matrix (Bi-nj) which contains 
neutrino mixing matrix is not known in details. 
Only t he upper limits, (s^')^ = T.j=4,- \Bun,\'^ 
[2012 1| . on the combinations of its elements are 
limited experimentally |22| . The leading-order 
LFV contributions come from the loops with two 
or four Bin matrices (there is one tree-level con- 
tribution, but it is negligible compared to the loop 
contributions). The combinations of matrix ele- 
ments which appear in the loop contributions to 
the LFV amplitudes cannot be determined, too. 
Therefore, for given (s^')^-s only the upper limits 
of the LFV amplitudes can be found. 

In model B the main sources of LFV are 
lowest-order loop amplitudes with two or four 
lepton-slepton-chargino/neutralino vertices |llj . 
The constants in these vertices are defined by 
mSUGRA SUSY breaking mediation mechanism 
|12| and the underlying GUT model, here cho- 
sen to be the minimal 50(10) model j4l5l6l7l8l^ . 
The model contains two SO{10) Higgs super- 
fields, 10 and 126 in the Yukawa sector, and 
additional Higgs superfields 126 and 210 to pre- 
serve the supersymmetry at the GUT scale 
and to break the 50(10) symmetry to 5C/(3)c x 
SU{2)l X f/(l) (G321) symmetry pi^ . By 
assumptions of the mSUGRA mediation mech- 
anism, the soft SUSY breaking terms at the 
GUT scale (the mass term and the trilinear in- 
teraction term) satisfy the universality condi- 
tions, according to which the mass soft SUSY- 



breaking term is flavour-diagonal. In the renor- 
malization group (RG) evolution from the GUT 
scale to the SM scale, mass soft SUSY-breaking 
terms receive the off-diagonal LFV contributions 
through the Dirac-neutrino Yukawa couplings 
|10j . Diagonalization of the sfermion and gaug- 
ino mass matrices leads to the chargino and neu- 
tralino mass eigenstates with LFV lepton-slepton- 
chargino/neutralino vertices. Only three param- 
eters, defining the universal soft SUSY-breaking 
parameters at the GUT scale, have to be intro- 
duced in this procedure. Concerning the mixing 
matrices in model B, which define the mass eigen- 
states, they are known in details. Therefore, for 
a given set of initial parameters, LFV amplitudes 
can be precisely evaluated. That allows to deter- 
mine both the lower and upper bounds of LFV 
decay rates j211- 

Concerning the RGE procedure there are two 
subtleties. First, as there are only two Higgs su- 
perfields in the Yukawa superpotential, the lepton 
and quark masses at the GUT scale are defined 
in terms of two mass matrices. That leads to two 
equations relating the quark and lepton masses 
at the GUT scale, which depend on two complex 
parameters The quark and lepton masses at 
the GUT scale are obtained RG-evoluting them 
from the experimental values at the scale of the 
Z-boson mass. Inserting them into the two equa- 
tions one can find the solutions for the two com- 
plex parameters (the phase of one parameter can- 
not be determined). With that, fermion mass 
matrices (and corresponding Yukawa matrices) 
are determined at the GUT scale, except for the 
right-neutrino mass matrix. To determine the 
right-handed neutrino matrix, one more complex 
parameter has to be introduced [5]. With right- 
neutrino masses determined, the light neutrino 
masses can be found using the type I see-saw 
mechanism. RG-evoluting the light neutrino op- 
erator one obtains the light-neutrino mass ma- 
trix at the low-energy scale which has to agree 
with the present neutrino-oscillation data. The 
agreement can be achieved. In fact, fitting of 
the neutrino-oscillation data is taken as a neces- 
sary condition that has to be fulfilled before pro- 
ceeding to further calculations. That means that 
in the LFV breaking couplings, the information 
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on the experimental neutrino oscillation data is 
built in. Second, the 50(10) model predicts the 
proton decay. The fitting of the model parame- 
ters can be achieved to satisfy the present pro- 
ton decay lower bound A detailed analysis 
of the Higgs mass matrices is needed to obtain 
the parameters defining the proton decay, and 
RG evolution of these parameters to the proton 
mass scale is necessary to obtain the the effective 
proton-decay Lagrangian |23I24| . Therefore, in 
principle, model B can relate branching ratios for 
neutrinoless LFV decays of charged leptons with 
the neutrino-oscillation data and lifetime of the 
proton. 

3. Effective Lagrangians for LFV interac- 
tions 

In any model containing SM as its low-energy 
limit, the following effective Lagrangians describe 
the neutrinoless LFV decays 

+a^,d-A^{Vt,PL+Vi^PB)%. (1) 
Af^.z = 9Ui,Z''{'PzPl+V^'Pr)%, (2) 

it^.H = gUrkPL + vSPR]m, (3) 

quQb XX' r 

x{ljTPxii){qaTPx'qb) (4) 

where Pl.r — (1 ± 75)/2. The dummy variables 
of the sums in the last expression assume fol- 
lowing values qaqb = uu,dd,ss,ds,sd, X,X' = 
LL, RR, LR, T = T^,!,^^!/- The Lagrangians 
contain the SM fields only. The effective La- 
grangians describe the decays £ ^'7*, £ £'Z* , 
£ £'H* and £ — > £'qaqb, respectively, where 
£ and £' , 7*, Z* and qa and qb denote leptons, 
photon, Z-boson, Higgs and quarks, respectively 
(star denotes an off- mass-shell particle). The ef- 
fective Lagrangian for the decay of lepton into 
three leptons is not given because we do not con- 
sider these decays. They have been considered 
previously in model A |18ll9j and B It has 
almost the same Lorentz structure as the effec- 
tive Lagrangian for £ £'qaqb- The informa- 



tion on LFV induced by a model induces is con- 
tained in the form factors 1^12^, Pi^2-yi Pz^ Pz^ 
Ph' Ph' ^Fab ■ Form factors contain the loop 
functions of loop diagrams describing the LFV 
processes. The structure of the "photon" effec- 
tive Lagrangian reflects the gauge invariance of 
the corresponding amplitude. The effective La- 
grangian for the "box" process £ — > £'qaqb is the 
most general one. In many models just few of 
the form factors are different from zero. Model A 
contains only terms bilinear in "ffj,PL,R matrices, 
while model B contains all form factors. 

Using the effective Lagrangians the neutrino- 
less SL LFV, amplitudes of charged leptons and 
amplitudes for £ — * ^'7 processes may be found 
at the quark-lepton level. 

4. Hadronization procedure 

To obtain the amplitudes of neutrinoless SL 
LFV charged lepton decays in terms of lepton 
and meson fields, the quark currents have to be 
transformed into the corresponding meson fields. 
Here only the decays with one or two pseudoscalar 
mesons or one vector meson are considered. The 
conversion of the axial-vector and vector cur- 
rents into pseudoscalar and vector mesons, re- 
spectively is achieved invoking the PCAC hypoth- 
esis [22] and VMD (vector-meson dominance) 
hypothesis j28l29l3()] . respectively. The scalar 
and pseudoscalar currents are hadronized com- 
paring the QCD quark Lagrangian with the cor- 
responding effective meson Lagrangian [21]. The 
tensor currents are hadronized using the equa- 
tions of motion for current quarks [321 ^^'^ VMD 
hypothesis. The semileptonic LFV amplitudes 
for processes with two pseudoscalar mesons in 
the final state have scalar-quark current con- 
tribution and vector-quark current contributions 
through the vector-meson resonances decaying 
into two pseudoscalar mesons. For that purpose 
the Lagrangian of vector-meson — pseudoscalar- 
meson interactions is introduced |33I34| . 

5. Comparison of the branching ratios 

With amplitudes at lepton-meson level, the 
comparison and evaluation of the decay rates 
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of the neutrinoless SL LFV charged lepton de- 
cays and i decays can be performed. For 
comparison of decay rates it is useful to clas- 
sify them in terms of the basic subamplitudes 
they contain, A^, Az, Ah and Abi^, which cor- 
respond to the effective Lagrangians (Q), |(5J, 
© and respectively. Such a classification is 
given in Table 1 for the processes considered here. 



Table 1 

The decay rates for the charged lepton neutri- 
noless SL LFV processes and I l''y processes. 



process 




model A 


model B 






7 


7 


(r ^ iP°] 




7. big 


7, Z, big 


(r ^ IP°] 


In 


l-q-box 


big 


(r -> 


\c 


7. Z, big 


7, Z, big 


(r ^ ^yo; 


In 


big 


big 




PlP2)c 


7, Z, H, big 


7, Z, H, big 


(j- I'' 


P1P2U 


big, W+W- 


big 


[t- l'- 


PlP2)c.H 


H, W+W- 


H 



c(n) denotes (non)conservation of the quark flavour 



In model A there is a window of the heavy 
neutrino masses (m^r) where the matrix elements 
BiM are not constrained by heavy neutrino mass 
(jun) values. Within that window, the functional 
dependence of the subamplitudes corresponds to 
the functional dependence of loop functions con- 
stituting the subamplitudes. SpecifHcally, the 
subamplitudes Ay, Az, Ah and Abi^ behave as 
Inmjv, iTij^, m?^ and InwAr, respectively. The 
Az has the largest combination of couplings and 
if it contributes to a LFV amplitude, it deter- 
mines its behaviour in the large mass limit. For 
very large masses, perturbative unitarity bound 
forces Biisi to behave as inverse of the heavy neu- 
trino mass leading to InTOjv/'Ti^, 
and In rnjv/^jv behaviour of the amplitudes A^, 
Az, Ah and Abi^ respectively. That behaviour 
assures that there is no LFV when the heavy neu- 
trino mass tends to infinity. Comparing the ex- 
pressions for the branching ratios in which only 
the dominant subamplitude is retained, one finds 
that there are only two groups of processes that 
are experimentally interesting, one with A^ only 
— > ^'7 processes), and the other with the 



Az amplitude (all quark-flavour conserving pro- 
cesses. The processes with box amplitude are 
strongly suppressed by their ttiat dependence and 
couplings in the vertices. The Higgs amplitude 
is strongly suppresses by the Higgs Yukawa cou- 
plings. Within the group with A^ amplitude, by 
far the best process is /x ^ 67. The branching ra- 
tios for the processes with the Az amplitude are 
all of the same order of magnitude, the t ^ tt", 
and T ^ p being somewhat better than the other 
processes. 

In model B, the characteristic scale of new 
physics is the mass of the lightest SUSY parti- 
cle TOsusY- Not only the loop functions in the 
Ay, Az and Abi have the softer behaviour than 
in model A, l/wgusY' I'^'^SUSY and l/Tilugy 
respectively, but all coupling constants decrease 
with msusY at least as fast as l/mgusY- That 
leads to the dominance of the Ay (photon penguin 
amplitude) in the large mass limit. Again, there is 
no LFV in the limit when the characteristic mass 
of the new physics tends to infinity. All branching 
ratios decrease as msuSY is increased, opposite to 
the behaviour of branching ratios in model A. We 
do not discuss the SM Higgs amplitude, because 
the model B has five Higgs bosons and some of 
them may have the LFV at the tree level. Limit 
on the tree-level Higgs contribution constraints 
combination of the Higgs boson mass and its cou- 
plings. In model B, the region of allowed values 
of mass parameter mguSY is small. As explained 
above, the couplings are precisely known when 
the parameters of the model are given. Therefore, 
upper and lower bound for tosusy can be deter- 
mined if there are adequate experimental data. 
The lower msuSY bound (tosusy — 560 GeV) is 
obtained from the /i — > 67 process. The upper 
bound msuSY — 800 GeV is determined from the 
present discrepancy between values of the muon 
magnetic moment obtained from the e^e~ data 
and T data |37j which cannot be explained within 
the SM at the moment. Within so small range of 
Tosusy values, it is hard to determine which of the 
amplitudes dominates without detailed numeri- 
cal study. The numerical values show that the 
different processes have different behaviour, with 
softer Tosusy for £ — > £'7 processes. That indi- 
cates that the contributions of some of three con- 
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tributions compete for dominance in the branch- 
ing ratio. Again, among all the processes the 
/i ^ 67 is the most interesting from the exper- 
imental point of view. All numerical results for 
both model A and B are shown in Table 2. For 
model B we find that in the parameter region, 
560 GeV < msuSY < 800 GeV consistent with 
all the experimental data, the predicted values for 
the semileptonic LFV branching ratios are close 
to the corresponding current experimental upper 
bounds. Therefore, they may be tested in the 
near future. 
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Table 2 

The branching ratios for neutrinoless SL LFV of a charged lepton and and t £'j processes. 



decay 


modes 


exp. bound 


model A 


model B 








'1 


1.2 X 10"" 


8.1 X 10-9^2^ 


1.2 X 10- 


11 




T ^ 


'1 


3.8 X 10-^ 


3.4 X IQ-^xl^ 


(2.6 X 10 


-10 


1.0 X 10-^) 


T ^ fl- 




3.1 X 10-'^ 


6.7 X 10-^4^ 


(l.l X 10 


-9 


4.5 X 10-^) 


T ^ e" 




1.9 X 10-'^ 


2.8 X 10-6y2^ 


(7.4 X 10 


-9 


4.1 X 10-*) 


T ^ fl- 




4.3 X 10-^ 


5.4 X 10-7y2^ 


(3.4 X 10 


-8 


1.8 X lO-'^) 


T ^ e" 




2.3 X 10-^ 


4.0 X lO-'^y^^ 


(l.l X 10 


-9 


6.2 X 10-^) 


T fl 




1.3 X 10^^ 


7.8 X 10-«y2^ 


(5.1 X 10 


-9 


1.3 X 10-*) 


T 


-p" 


2.0 X 10"*^ 


2.7 X 10-6y2^ 


(4.0 X 10 


-14 


2.9 X lO-i'^) 


r ^ 


-p° 


6.3 X 10-*^ 


5.3 X 10-^2/2^ 


(1.3 X 10 


-13 


1.2 X 10-12) 


T ^ e' 


^0 


6.9 X 10^6 


2.7 X 10-6y2^ 


(9.6 X 10 


-10 


5.1 X 10-^) 


T fr 


-0 


7.0 X 10-6 


5.3 X 10- 


(4.1 X 10 


-9 


2.2 X 10-*) 


T 




2.2 X 10-6 


2.7 X 10-6y2^ 








T ^ fl- 




8.2 X 10-6 


5.2 X 10-7y2^ 








T ^ e" 


K+K- 


6.0 X 10-6 


1.1 X IQ-^yl^ 








T fj,- 


-K+K- 


1.5 X 10-5 


2.1 X 10-^^ 









The model A values are upper limits of branching ratios multiplied by combinations of Bin matrix 
elements that are always smaller than one [3^1 ■ For model B the first (second) value in brackets represents 
lower(upper) bound for a branching ratio. For /i — > 67 only the upper bound is given, which served as 
an input parameter to determine the minimal msusY value. 
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